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Abstract. Using the USNO CCD Astrograph all-sky Catalog (UCAC2), we measure the mean proper motion of the two 
Magellanic Clouds. Appropriately-selected LMC populations show a proper motion (p.a,p.S) - (+0.84, +4.32) that is signif- 
icantly higher, in {/id), than currently accepted Hipparcos-like values; {fj.a,fj.S) - (+1.94, -0.14). A higher (p.a) value is also 
found for the SMC. Interestingly, the mean UCAC2 LMC proper motion agrees very well with the only work in the literature 
(Anguita et al. 2000 1 pointing to an unbound Magellanic Clouds-Milky Way interaction. Nonetheless, the implications of the 
UCAC2 proper motion are hard to reconcile with our present day understanding of the Clouds-Galaxy interaction unless one 
assumes a more massive Milky Way. Consequently, although no sources of systematic error have been identified, it is perhaps 
most likely that the UCAC2 catalog has an as yet unidentified systematic error resulting in an inconsistency between UCAC2 
and Hipparcos based results for the Magellanic Clouds. 
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1. Introduction 

The gas envelope surrounding both the Large and Small 
Magellanic Clouds (LMC and SMC) is reminiscent of past in- 
teractions between the two galaxies. In turn, the Magellanic 
Stream is evidence of past interactions of both galaxies with 
the Milky Way (Kroupa & Bastian [l997t . Modeling the past 
stellar and dynamical evolution of the Magellanic Clouds re- 
quires good knowledge of many parameters, in particular those 
concerning their spatial velocity vector and transverse veloc- 
ity (v,, van der Marel 2004). A sound determination of the lat- 
ter parameter holds the key for a better understanding of the 
Milky Way-Magellanic Clouds interaction. The LMC/SMC 
transverse velocity can be measured directly from the observed 
proper motion of the galaxies, or alternatively, via (i) dynam- 
ical modeling of the Magellanic Stream or (ii) analysis of the 
line-of-sight velocity field of the LMC/SMC. 

Many groups have addressed the v, measure of the 
Magellanic Clouds via the direct method of proper-motion 
(pm). Being closer, the LMC has been the subject of more de- 
tailed studies (see Table ^ f° r a compilation of pm values). 
At first sight, one notes an overall excellent agreement be- 
tween the first 6 determinations. However, there is one outlier 
(~ 5.5<t times the average value in declination), namely that 
of Anguita, Loyola and Pedreros (2000 hereafter ALP00). The 
ALP00 study is based on 125 CCD images spanning 8 years, 



and put on a QSO absolute reference frame. It was the first pm 
analysis of LMC stars reaching V - 22; i.e. including nearly 
all its composite stellar populations. 

In a later work, Pedreros et al. (2002 1 using a similar CCD 
data set and the same QSO method, contradicted their earlier 
ALP00 results and derived Hipparcos-cons/sfenf pm values. 



Puzzled by this discrepancy, Pedreros et al. yP02i searched 
for systematic errors possibly affecting the ALP00 study. Their 
most plausible explanation pointed to peculiar pm or peculiar- 
ities in the rotation curve of the LMC between their two stud- 
ied fields (located at different positions and distances from the 
LMC center). Further progress on this issue awaited results 
from the 9 million LMC stars of the MACHO project. To our 
knowledge, the only MACHO study on the LMC pm is that of 
Drake et al. (|2D02j, confirming "normal" pm values. Thus, the 
ALP00 study remains the only, inconsistent, pm value pointing 
to the Magellanic Clouds being unbound. 
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Following the methodology outlined in Momany et al. |2004 1 
we analyze UCAC2 (Zacharias et al. 2004 1 pm data of the 
Magellanic Clouds. Interestingly, we find that the UCAC2 data 
unequivocally point toward a high (fid) of the LMC, as in the 
ALP00 study. In this paper we attempt to assess the reliabil- 
ity of the UCAC2 pm data by searching for color/magnitude 
dependencies or other systematic effects. The UCAC2 pm 
data show no significant evidence of systematic errors. The 
Hipparcos-UCAC2 discrepancy basically remains unsolved. 
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Table 1. A compilation of LMC proper motion studies. 



Source 








System 




(mas yr ) 


(mas yr~') 




Kroupa et al. (1994) 


+ 1.30:1 


:0.60 


+ 1.10±0.70 


PPM 


Jones etal. (1994) 


+ 1.20:1 


:0.28 


+0.26±0.27 


Galaxies 


Jones etal. (1994) 


+1.374 


:0.28 


-0.18±0.27 


Galaxies 


Kroupa Bastian (1997) f 


+ 1.94±0.29 


-0.14±0.36 


Hipparcos 


Drake et al. (2002)+ 


+ 1.40±0.40 


+0.38+0.25 


MACHO 


Pedreros et al. (2002) 1 ' 


+2.00:1 


:0.20 


+0.40±0.20 


QSO 


Anguita et al. (2000) 


+ 1.70:1 


:0.20 


+2.90±0.20 


QSO 



1 The reported pm are those corrected and listed in van der 
Marel (20021. 



The plausibility and implications of the UCAC2 and ALP00 
values is discussed. 



2. The UCAC2 LMC and SMC proper motion 

The recent availability of large, high precision astrometric and 
photometric catalogs like 2MASS (Cutri et al. 120031 . SPM3 
(Girard et al. l2"U04T l and UCAC2 (Zacharias et al. l2"U04l al- 
low kinematic investigations of structures within the Galaxy 
and determination of the orbits of its nearby satellite compan- 



ions. In Momany et al. (2004), using the UCAC2 data, we 
have been able to recover/confirm the HST pm measurement 
of the Sagittarius dwarf galaxy located at a distance of 25 kpc: 
(pa,p$) ^ (-2.5,-2.0) mas/yr (Irwin et al. 1996). This value 
has also been recently confirmed by Dinescu et. al. d2004> us- 
ing the SPM3 catalog. Thus, the Sagittarius dwarf is an ex- 
cellent example of the reliability of UCAC2 pm data, and en- 
courages the use of these data for a similar analysis for the 
Magellanic Clouds. 

Combining the 2MASS infrared photometry with UCAC2 
optical and pm data offers a unique opportunity to discriminate 
LMC/SMC stellar populations from the Galactic (disk/halo) 
contribution. On the one hand, this will help identify possi- 
ble dependencies of the pm data on the employed photometric 
magnitudes and colors, and on the other, it dramatically ex- 
pands the number of pm stars in different LMC/SMC popula- 
tions (having different colors and magnitudes). This reduces the 
error budget and possible contamination effects. The net result 
is a more robust estimate of the bulk pm. Zacharias et al. esti- 
mate UCAC2 pm errors to be 1 - 3 mas/yr for 12'' 1 magnitude 
stars and about 4-7 mas/yr at 16 th magnitude. Throughout this 
paper, we use absolute pm data for LMC and SMC stellar pop- 
ulations having 12 < R < 16. Following a 3cr clipping, we esti- 
mate the mean rms of the average pm values to be < 4.0 mas/yr. 
Although quite significant, such an rms error is acceptable con- 
sidering the distance of the Magellanic Clouds. Moreover, we 
prefer to assume such a high error rather than that obtained by 
dividing it by the square root of the number of stars. 

The first epoch data for the pm come from different cat- 
alogs including Hipparcos/Tycho, AC2000.2, as well as re- 
measured AGK2, NPM and SPM plates. As in Momany et 



al. (120041 ) . we limit the UCAC2 pm data to only one first epoch 
catalog in order to limit internal inhomogeneities, namely the 
AC2000.2 catalog (Urban et al. [T9981 l200Tt . Based on the 
Carte du Ciel Astrographic plates and measurements, the con- 
struction of the AC2000.2 was performed scanning all the col- 
lected plates and then collecting all the positional measure- 
ments in a homogeneous system based on Hipparcos. A key 
feature of AC2000.2 is the mean epoch of the positions which 
range from 1890 to 1940. In the particular case of the LMC the 
mean epoch of AC2000.2 is =s 1893.5. Comparing the LMC 
mean epoch of AC2000.2 with the observational epoch of the 
UCAC2 (1998. 1+2002. 9) gives a very large epoch span for the 
pm measurements, therefore reducing the pm errors. 

Figurensummarizes our analysis in deriving the LMC (up- 
per panels) and SMC (lower panels) mean pm. The left panels 
of Fig. □ show combined UCAC2/2MASS R, (R - K) color- 
magnitude diagrams, for two fields (1° radius each) centered 
on the LMC and SMC. Dark symbols are stars whose pm orig- 
inates only from the AC2000.2 catalog (and used in estimating 
the pm), whereas light symbols have pm from a number of dif- 
ferent catalogs, including Hipparcos. The middle panels high- 
light the selected LMC/SMC stellar populations. Relying on 



the color-magnitude diagrams in Nikoleav & Weinberg (2000) 
and Alcock et al. (2000) of the LMC, one can easily disentan- 
gle LMC stellar populations from the foreground Galactic con- 
tribution. Indeed, the two left-most selected populations (open 
triangles) are almost certainly LMC young main sequence and 
Helium burning stars. On the redder side, the three right-most 
populations (open circles) are red super-giants, oxygen- and 
carbon-rich LMC stars, respectively. In between the LMC pop- 
ulations, one identifies two vertical sequences (open squares) 
which are Galactic disk F-K dwarfs and K giants. As noted by 
Alcock et al. (2000), in these regions of the color-magnitude 
diagram the Galactic disk stars out-number those belonging to 
both the Thick Disk and Halo combined. A confirmation of this 
comes from a Besangon synthetic color-magnitude diagram of 
the Galactic contribution: modeling in both the LMC and SMC 
directions clearly showed that virtually no Galactic stars fall 
in the Magellanic Clouds selected regions. Filled circles are 35 
Hipparcos stars tabled in Kroupa & Bastian ( 1997 their Table 5 
excluding star number 26222). Hipparcos stars plotted as open 
circles are those falling inside our 1° radius UCAC2 fields. 

The right panels show the averaged 3<x clipped (pa, pd) di- 
agram of all selected populations. For the LMC (upper panel) 
we easily distinguish three separate groups, namely: (i) the 
Galactic populations characterized by a high (pd) - 8.5 
mas/yr; (ii) the LMC populations clearly clustering around 
(pa,pd) ^ (0.9,3.9) mas/yr; and (Hi) the Hipparcos mea- 
surements at (pa,p5) - (1.5,-0.3) mas/yr. Likewise for the 
SMC (lower panel): (i) the Galactic populations at (pa,p6) ^ 
(9.0,-1.5) mas/yr; (ii) the SMC populations at (pa,p5) - 
(4.5,-3.3) mas/yr; and (Hi) the Hipparcos measurements at 
(pa,p6) (1.0, -1.0) mas/yr. The right panels of Fig^show 
the discrepancy between the Hipparcos and UCAC2 pm mean 
values, mainly in (pd) for the LMC and in both directions for 
the SMC. Overall, the UCAC2 data indicate a mean pm for 
the LMC of (pa,p6) - (0.84, +4.32) mas/yr based on ~ 1620 
stars and ^ (4.44, -2.92) mas/yr based on ~ 1200 stars for the 
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SMC. Lastly, we searched for possible variations of the mean 
pm across the LMC. However, for all considered fields, the re- 
sulting pm in declination was always larger than Hipparcos- 
like values. In particular, we were able to recover/confirm the 
ALPOO pm value of their 2 fields centered at 5 h 57"' and 5 h 58'". 
Thus, the inconsistency of UCAC2 data is found all over the 
LMC. 

One may be surprised by the high pm values obtained for 
the Galactic populations. In this regard, we note that the de- 
rived pm of Galactic populations is very well reproduced by 
the Besangon Galactic model simulations. The mean model 
pm was computed on the synthetic color-magnitude diagrams 
for the two vertical sequences and is plotted as an open-starred 
symbol in both the right-most panels. We emphasize that these 
are model expectations, meant only to test the reliability of 
UCAC2 values for the disk stars. The agreement between 
UCAC2 and Besancon however remains quite significant. 

3. The UCAC2-Hipparcos discrepancy 

Unveiling the origin of the UCAC2-Hipparcos discrepancy is 
quite a difficult task since it may be due to many different 
sources of systematic errors. In the following we list a num- 
ber of possible causes. 

- Magnitude/color dependencies: The LMC/SMC 
UCAC2/2MASS selected populations span ~ 5.0 magni- 
tudes in (R - K) showing no appreciable dependence on 
color. Furthermore, dividing the 5 selected populations in 
magnitude bins did not result in any appreciable variation 
of the mean pm. Hence, although the presence of residual 
color/magnitude terms in any pm analysis is almost 
unavoidable (Platais et al. 2003 1, these seem not to be 
significant enough to explain the discrepancy; 

- Zone systematics: In Zacharias et al. (2004) pm compar- 
isons in the LMC field were computed for a sample of 
each of 200 stars. This star by star comparison yielded a 
difference of (Apa,Ap5) = (+0.2,-1.0) mas when com- 
pared to the average value of van der Marel et al. (|2002i. 
This is clearly not enough to explain the (Apd) 4.0 mas 
Hipparcos-UCAC2 difference. Zacharias et al. also checked 
for systematic differences in the UCAC2 data as a function 
of declination zone (their Table 8), again not showing any 
significant systematic effect; 

- Galactic contamination in the Hipparcos sample: The 
Hipparcos sample has been selected on the the basis of 
radial velocity. These are bright, young O/B LMC/SMC 
members which on Hipparcos color-magnitude diagrams 
(see also Fig. 1) are clearly limited to the brightest regions 
(R < 12.00), and in a narrow color range < (R — K) < 
1.5. In this color-magnitude region some contamination 
by nearby (within 1 kpc) disk main sequence stars is pre- 
dicted in the Besangon photometric/kinematic simulation. 
However the expected radial velocity of these simulated 
stars is well below 100 km/s implying that no contaminants 
should be present in the Hipparcos LMC/SMC sample; 

- Hipparcos O/B LMC stars selection effects: The Hipparcos 
LMC stars (V > 10.5) are actually far below the 100% com- 



pleteness magnitude limit of Hipparcos (V - 9), and are 
located in a color-magnitude range where the completeness 
rapidly falls to zero around V - 13. Moreover, these stars 
are projected on the main body of the LMC and are sur- 
rounded by a rich and crowded background that could lead 
to an incorrect centroiding measurement; 
- O/B spectral type dependencies in Hipparcos: A recent de- 
tailed analysis by Schroeder et al. (2004) of the Hipparcos 
parallaxes of Galactic O/B stars pointed out that these are 
heavily erroneous leading to absolute magnitude estimates 
up to 5 magnitudes fainter. On the other hand, Pan et al. 
(2004) identified a duplicity in the star Atlas affecting the 
centroiding of Hipparcos, and suggested this as the cause 
of the anomalous (10% smaller than accepted) Hipparcos 
distance for the Pleiades. A bias in the parallaxes of bright 
O/B stars in the Pleiades has also been confirmed recently 
by Soderblom et al. (2004). It is not clear to us how wrong 
parallaxes can affect proper motion determination, however 
we note that Hipparcos gives non-negative parallaxes for 
most of the LMC stars. This is much larger than the ex- 
pected value for the LMC. 

In summary, we are unable to find a consistent explanation 
as to why there exists a discrepancy between the Hipparcos and 
UCAC2 pm measurements. Any a priori comparison of the two 
catalogs would certainly favor the UCAC2 as the data source 
on which to base a pm determination of the LMC. Indeed: (i) 
the large number statistics involved in the UCAC2 LMC/SMC 
measurements; (ii) the large time baseline for the pm measure- 
ments; (Hi) the reliability of UCAC2 in determining the pm for 
the Sagittarius dwarf; ( iv) the excellent consistency of the mean 
pm of different LMC/SMC populations with different colors 
and magnitudes; and (v) the good reproduction of the mean pm 
of the Galactic populations in the LMC/SMC fields, are all evi- 
dence in favor of the UCAC2 measurements with respect to the 
limited (in magnitude and spectral class) Hipparcos sample. 

4. Discussion: implications of the UCAC2 proper 
motion 

Serious problems arise when one replaces the Hipparcos with 
the UCAC2 value for the LMC proper motion. Indeed, the 
physical implications of UCAC2/ALP00 pm contrast with 
our present-day understanding of the Local Group in gen- 
eral, and the Cloud-Galaxy interaction in particular. All meth- 
ods used to measure the Clouds' transverse velocity (v,) ar- 
gue against the UCAC2 pm values. Applying eq. (55) of 
van der Marel et al. (2002) for a Hipparcos-like pm value, 
the resulting three-dimensional space velocity for the LMC 
is I^lmc = (-56,-219,186) km/s, while the corresponding 
total, tangential and radial velocities are (293,281,84) km/s, 
in the Galactocentric rest frame. On the other hand, apply- 
ing the UCAC2 LMC pm gives "v^lmc = (-974, -195, -113) 
km/s, with the same velocities being (1001,971,241) km/s. 
Similarly, the UCAC2 pm of the SMC would imply ~^smc = 
(-634,-974,475) km/s, and (1256,1251,103) km/s. These 
UCAC2 -based values are hard to reconcile with our current un- 
derstanding of the: 
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- LMC rotation: the UCAC2 implied tangential velocity of 
the LMC, 971 km/s, is so large that, if present, should ap- 
pear as in a prominent solid-body rotation component as 
expected from the large projected dimensions in the sky of 
the LMC. There is no evidence of such a high component 
(see van der Marel 2002); 

- Local Group: in the Galactocentric rest frame, there is 
no evidence in Local Group or nearby groups of galax- 
ies of such extreme velocities. As recently shown by 
Karachentsev (2005 1, the velocity dispersion of dwarf 
galaxies under the influence of the Milky Way is cr vr st 
86 km/s. Similar values of the velocity dispersion are found 
in other nearby groups like Andromeda, CenA, M81/82 
etc., indicating that no peculiar motions are present. The 
LMC radial velocity based on UCAC2 is almost 3x this 
typical velocity dispersion, i.e. much higher than the escape 
velocity; 

- Magellanic Stream: the UCAC2 LMC transverse veloc- 
ity puts it at ~ 90° from the direction of the Magellanic 
Stream. Although we still do not know why the stream is 
starless, it is hard to envisage how this relatively young fea- 
ture survives at right angles to the motion of the LMC (see 
Mastropietro et al. 2004 for recent modeling); 

- Magellanic Bridge: analysis of the Magellanic Stream and 
Bridge shows evidence that the Bridge is old (Putman et al. 
2003 1. This strengthens the idea that the LMC and SMC 
have jointly orbitted around the Galaxy at least for the past 
Gyr or so. However, this observational evidence is hard 
to reconcile with the hypothetical UCAC2 scenario where 
there exists a relative velocity between the LMC and SMC 
of the order of 1033 km/s. 

A scenario with the LMC traveling at a high velocity (first 
passage) in the vicinity of the Milky Way seems also to con- 
tradict recent observational (Cole et al. 2005 1 and theoretical 
(Bekki & Chiba 2005 1 studies that successfully reproduce the 
observed morphology, dynamics and formation histories of the 
LMC assuming it has been in tidal interaction with the SMC 
and the Milky Way. The only viable scenario that allows for 
the UCAC2/ALP00 LMC pm is one in which the Milky Way is 
more massive (within 50kpc) than 3 H-4x 1O 12 M as suggested 
by Mendez et al. {T999 l. Interestingly, the Mendez et al. mod- 
eling is based on a proper motion analysis of 30,000 stars from 
the Southern Proper-Motion (SPM) survey, and is consistent 
with the Anguita et al. pm value of the LMC. Nevertheless, the 
Mendez et al. Milky Way total mass estimate remains much 
larger than those reported in recent studies, e.g. Sakamoto, 
Chiba & Beers (2003 1 who estimate the mass of the Galaxy 
(within 50 kpc) to be 5.5 x 1O U M . 

In conclusion, the UCAC2 data may have been valuable in 
determining the pm of the Sagittarius dwarf at 25 kpc and in 
the direction of the Canis Major over-density at ^ 8 kpc, yet 
the UCAC2 pm data are hard to reconcile with our present day 
understanding of the Cloud-Galaxy interaction. Given that our 
study has failed to identify sources of systematic errors, the 
UCAC2/Hipparcos inconsistency remains unsolved and awaits 
high-accuracy HST, GAIA and SIM proper motion measure- 
ments. 
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Fig. 1. Le/f panels display UCAC2/2MASS R, (R - K) color-magnitude diagrams of the LMC and SMC. Dark symbols are used 
for stars whose pm originates only from the AC2000.2 catalog, while light symbols have a pm coming from a number of catalogs, 
including Hipparcos. Middle panels highlight the extracted populations: the two vertical sequences around (R-K) a 1.5 are 
Galactic disk stars, the rest are LMC/SMC populations. Dots brighter than R ~ 12.0 are Hipparcos samples and tabled in Kroupa 
& Bastian ( 1997 1. Those plotted as open circles fall within the examined 1° radius UCAC2 field. Right panels show the mean 
proper motion of the selected populations in the (pa,p6) plane. The different symbols correspond to those in the middle panels. 
To these we added, an upside down open triangle marking the ALP00 value, while an open starred symbols shows the expected 
mean pm of Galactic disk stars (as derived from the Besancon simulation). 



